Possibility of "magic" co-trapping of two atomic species in optical lattices 
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Much effort has been devoted to removing differential Stark shifts for atoms trapped in specially 
tailored "magic" optical lattices, but thus far work has focused on a single trapped atomic species. 
In this work, we extend these ideas to include two atomic species sharing the same optical lattice. 
We show qualitatively that, in particular, scalar J = divalent atoms paired with non-scalar state 
atoms have the necessary characteristics to achieve such Stark shift cancellation. We then present 
numerical results on "magic" trapping conditions for ^^Al paired with *^Sr, as well as several other 
divalent atoms. 
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I. INTRODUCTION 

Cold atoms and molecules have proven to be a use- 
ful platform for a variety of precision measurement and 
quantum information experiments. Optical lattices pro- 
vide a method of trapping cold atoms and molecules that 
minimizes Doppler shifts and allows long measurement 
times. However, the trapping field causes shifts in the 
internal energy levels via the ac Stark effect, and these 
shifts depend on the intensity of the trapping lasers. This 
problem has been considered in detail in the context of 
optical lattice clocks (see, e.g., [il-|j|)- The solution em- 
ployed in these works is to use an optical lattice operating 
at "magic" trapping conditions (e.g., magic wavelength) 
for which the energy levels of interest experience the same 
shift, even as the atom or molecule moves within the lat- 
tice. This leaves the transition frequency of interest un- 
changed, as long as higher-order effects are negligible. In 
another recent work, the problem of removing Stark and 
Zeeman shifts simultaneously was considered for alkali- 
metals |5j- In each of these works, cancellation of the 
Stark shift was achieved for a specific clock or qubit tran- 
sition in a single atomic species. 

In other experiments, simultaneous trapping of two 
atomic species has been achieved using double magneto- 
optical traps (MOTs) [aTgl, combined with an optical 
dipole trap in 8] , with the primary goal of sympatheti- 
cally cooling atomic species that are difficult to cool di- 
rectly. However, in these experiments the internal states 
of the atoms are heavily perturbed by the MOT, which is 
disadvantageous for precision measurements. Comagne- 
tometry experiments also utilize two atomic species shar- 
ing the same volume. One such apparatus was used to 
create a highly accurate nuclear spin gyroscope [9| , while 

^ We use atomic units except where otherwise noted. 



other groups have used comagnetometers in searches for 
physics beyond the standard model, such as setting up- 
per bounds on P,T-violating permanent electric dipole 
momei it i lO.1 as well as bounds on Lorentz and CPT vio- 
lation [ll|~^ 

In the present work, we aim to merge these ideas 
by considering two atomic species trapped in the same 
"magic" optical lattice. Our goal is to achieve Stark- 
shift cancellation simultaneously in both atomic species, 
for a single two-level transition in each. We first intro- 
duce the equations underlying the ac Stark effect. These 
equations guide our selection of favorable atom pairs to 
analyze in detail. Finally, we present results of numer- 
ical calculations on ^^Al paired with Sr, Yb, and other 
divalent atoms. 



II. STARK SHIFT 

The derivation of the Stark shift follows closely the 
treatment in [3| and [12| , so only a summary of the results 
will be given here. The ac Stark shift of an atomic state 
\F, Mp) with total angular momentum F and projection 
Mp can be written as ^ 
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where £l and lo are the amplitude and frequency of 
the lattice light field. The frequency-dependent atomic- 
structure pre-factor a^§\,j (w) is the total ac polarizabil- 
ity of the state; it may be decomposed in terms of irre- 
ducible tensor operators as 



J 



tot 
'-'FMf 



(C) 



aUuj) 



{k 



,)A^aUoj) 



l{ 



Sle-e. 



3Mf ^ - F{F + 1) 

F{2F-1) 



T 

ap 



(c) 



(2) 



Here A is the degree of circular polarization, af., a^, 
and a^ are the scalar, vector, and tensor polarizabilities, 
respectively, and the unit vectors are the lattice laser 
wavevector (k), polarization (e), and quantization axis 
(e^). Fig. [1] depicts the geometry. The quantization axis 




is defined by an applied static magnetic field, whose ori- 
entation relative to the optical lattice may be arbitrary. 
This "quantizing magnetic field" ensures that Mp re- 
mains a "good" quantum number for the ac Stark-effect 
perturbation formalism, from which Eq. ([2]) is derived; 
the energy shifts due to the Stark effect must be small 
compared with the Zeeman splitting of the magnetic sub- 
levels. It is important to note that the total polarizability 
depends not only on the spatial orientation of k, e, and 
Cz, but also on the lattice frequency uj as well. Also use- 
ful is the total differential polarizability Aa between two 
states \F',M'p) and \F,Mp), given simply by 
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Analogously, we will use the scalar differential polariz- 
ability, defined as 



Aa^{uj) = a%,M, iuj) ~ af-^A^) 



(4) 



FIG. 1: (Color online) Relation of angles to unit vectors for 
linearly polarized light, k, e, and e^ are the lattice wavevector, 
lattice polarization, and quantization axis, respectively. For 
circular polarization, we take i to be complex, following the 
Jones calculus conventions. In this case dk is unchanged, while 
6p is no longer well-defined. 



and also the vector and tensor differential polarizabilities, 
defined similarly. 

For linearly polarized light, the vector contribution 
vanishes (^ = 0), and e-e^ = cos Op, where 9p is the angle 
between the polarization and quantization unit vectors. 
Then the total polarizability becomes 
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In the case of circularly polarized light, k ■ iz — cos 9k 
and |e-ez| — 5 sin 9k, where Ok is the angle between 
the wavevector and quantization unit vectors, e is taken 
to be complex and constant in time, rather than real 



and time varying. This means that, although counter- 
intuitive, le-e^l is time-independent even for circular 
polarization. The total polarizability is then 
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It should be noted that in both cases, the Stark shift 
depends on only one angle (either 9p or 9k , but not both) . 
This means that for purely circularly or linearly polarized 
light, we have only two "handles," not three, with which 
to achieve "magic" conditions: a single angle and the 
lattice frequency uj. 

It should be noted that the tabulated formulae are 
valid only for a sufficiently weak magnetic field. For finite 
B fields, the Stark shift acquires an additional cross-term 



formed as a product of the Zeeman and vector light shifts; 
this provides an additional knob for tuning magic condi- 
tions (see [ 13l. Il4l| ) : we will comment on this in SectionFvl 



III. SELECTING ATOMS FOR TRAPPING 

Guided by the above equations, we state some general 
juidelines in selecting two atomic species to trap. Con- 



sider first the vector and tensor polarizabilities. These 
are rank 1 and rank 2 irreducible tensor operators, re- 
spectively. Triangular selection rules dictate that a|^ 
must vanish for F — 1/2 states, and that both a^ and 
ap must vanish for F — states. Alternately, the po- 
larizabilities of F = states have no dependence on the 
angles 0p or 9k- 

This leads us to consider divalent atoms, such as Sr or 
Yb, which have ^^o ground states. The bosonic isotopes 
have nuclear spin / = 0, so that F — J — 0: due to 
the angular selection rules a§^ and ap must vanish. The 
fermionic isotopes have / ^ 0, and so (for / > 1/2), a]^ 
and ap do not vanish identically. However, a-^ and ap 
are nonzero only due to a highly suppressed contribution 
from the hyperfine interaction (HFI), and so ap domi- 
nates. For optical transitions between J = electronic 
states, the differential polarizability is almost entirely due 
to the scalar polarizability. Therefore the light shift, be- 
ing of scalar nature, does not depend on atomic spatial 
orientation. We conclude that for — clock transitions 
in divalent atoms, the dependence on geometry is negligi- 
ble and "magic" conditions are set entirely by the lattice 
frequency. If one of our two trapped atomic species is a 
scalar atom, we must use a "magic" lattice frequency for 
that atom. Since the frequency becomes fixed, we have 
to achieve magic conditions for the second species solely 
by manipulating one angle, 0p or 0k (for fixed circular or 
linear polarization). 

While the discussion so far has been of a general na- 
ture, for illustration let us consider microwave hyperfine 
transitions attached to a J = 1/2 electronic state in our 
second atomic species. The relevant atoms may include 
the ground state alkali-metal atoms (si/2 ground states) 
or group III atoms {pi/2 ground states). To minimize 
Zeeman shifts, we will work with Mp = states exclu- 
sively, which also eliminates the vector term in Eq. ([B]). 
There are some technical but important peculiarities as- 
sociated with computing Stark shifts between levels of 
the same hyperfine manifold: we are considering the 
Stark shift of hyperfine levels attached to the same elec- 
tronic state. To the leading order, the shift is determined 
by the properties of the underlying electronic state. How- 
ever, because the electronic state for both hyperfine lev- 
els is the same, the scalar Stark shift of both levels is 
the same. An apparent difference between the two clock 
levels is caused by the hyperfine interaction (HFI), and 
the rigorous analysis involves so-called HFI-mediated po- 
larizabilities Q. Similarly, when neglecting the HFI, the 
tensor shift for the J = 1/2 state vanishes due to the an- 
gular selection rules; we are also led to a necessity to treat 
the tensor shift using the HFI-mediated polarizabilities. 

As demonstrated in Ref . [4| , for Mp — clock or qubit 
states and fixed B||A; or B _L fc geometry neither Rb nor 
Cs can be magically trapped. However, the magic trap- 
ping can be achieved for Al and Ga. Our consideration 
here is more general, as we also consider varying rotation 
angles. Still even in our more general setup we do not 
find magic angles for the Mp = clock or qubit states 



in Cs and Rb (again this statement is limited only to the 
weak B field regime). We focus below on Al, where the 
required magic angles can be found. 

To summarize, we choose a scalar divalent atom and 
a non-scalar valence state atom to trap simultaneously. 
We choose for our optical lattice a "magic" frequency for 
the divalent atom, so the differential polarizability for its 
optical transition of interest is zero. At this lattice fre- 
quency, the differential polarizability for the microwave 
transition in the non-scalar atom is non- vanishing, in gen- 
eral. But by tuning 0p or 0k (depending on lattice po- 
larization), we aim to create "magic" conditions in the 
non-scalar atom as well. 

We perform numerical calculations to find "magic" 
conditions for particular pairs of atoms. Utilizing a va- 
riety of relativistic many-body techniques, as described 
in [1^ , we calculate the atomic polarizabilities in Eq. ^ 
and Eq. © . We then numerically search for values of 0k 
or 9p that will make Aa vanish in the non-scalar atom. 
We focus our attention on ^''Sr paired with ^''Al, although 
we also present results for ^^Al paired with other diva- 
lent atoms. In ®^Sr, the transition of interest is between 
the 5s^ -^^o ground state and the 5s5p '^Pq metastable 
state [l^. In ^^Al, the transition of interest is between 
the \F = 2,Mp ^ 0) and \F = 3,Mp = 0) hyperfine 
states, both in the 3pi/2 electronic ground state. We 
work with Mp — states exclusively in Al so that (to 
first order) Zeeman shifts vanish. 



IV. RESULTS 

In Fig. [2] we plot results representative of our numeri- 
cal calculations. We graph the differential polarizability 
between the two hyperfine states in ^^Al as a function of 
angle, for linear and circular lattice polarizations. 

While section [mi suggests that Stark-shift cancellation 
may be possible for ^^Al at arbitrary lattice frequencies, 
it is not obvious if "magic" angles for ^^Al actually exist 
at the "magic" wavelengths for the divalent atoms of in- 
terest. Considering Eq. (O and Eq. ^, we note that the 
prefactor on the tensor term varies from approximately 
1 to — ^ for linear polarization and from 3: to — i for 
circular polarization. This means it is highly probable a 
"magic" angle exists if aj,, and a|^ are at least a few times 
greater than ap, and ap (barring a happenstance cancel- 
lation of a|^/ and aj,). Our calculations reveal that, over 
the range of wavelengths studied (see Table HI , aj^, and 
a|^ vary from 2 orders of magnitude larger than ap, and 
ap to roughly the same order of magnitude. Only near 
the Al resonance at 394 nm does the scalar polarizabil- 
ity dominate. This suggests that "magic" angles should 
exist for many pairings of Al with scalar divalent atoms. 

Table U shows results of our numerical calculations, 
confirming the above assertions. We tabulate the lat- 
tice frequency for which the Stark shift vanishes in each 
divalent species, and "magic" angles for which the Stark 
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We can write the fractional clock shift as 
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FIG. 2: (Color online) Differential polarizability (Eq. |3| be- 
tween \F = 3, Mf = 0) and \F = 2,Mf = 0) states in ^'^Al as 
a function of angle. The top graph is for a linearly polarized 
lattice, while the bottom graph is for a circularly polarized 
lattice. Note 9p is the angle between the polarization and 
quantization axis, while 9^ is the angle between the wavevec- 
tor and quantization axis, as shown in Fig. [l] The magic 
angle is the value of 9 for which Aa = 0, highlighted with 
red triangles. This is for a lattice at 813.42 nm, the "magic" 
wavelength for the optical clock transition in *^Sr. 



shift in ^^Al also vanishes for linearly and circularly po- 
larized lattices. The angles Op and Ok are as depicted in 
Fig. [T] and defined in Eqs. ([5]) and (H]), respectively. 



where II is the laser intensity and Aa denotes the differ- 
ential polarizability between the upper and lower clock 
levels, as defined in Eq. ^. For Al in a lattice with 
Sr (A = 813 nm), a misalignment of 0.1° corresponds to 
Aa w 10~^ a.u. If we assume a lattice intensity on the 
order of 10 kW/cm^, we are lead to a fractional clock 
shift of the microwave transition — « 10~^^. This effect 
would be somewhat worse for shorter wavelength lattices, 
since the polarizabilities tend to increase faster than v. 



V. RUBIDIUM AND CESIUM 

So far we focused our attention on the weak-field 
regime to reduce Zeeman-sensitivity of the M'p — Mp = 
transition. If we relax the weak-field requirement, magic 
conditions for such transitions can be found by varying 
the strength of the B-field. Magic trapping conditions 
then exist for all trapping wavelengths (barring atomic 
resonances), see Fig. 2 of Ref. [ij]. The required B field 
strengths are of a few Gauss. For example, one could read 
off Fig. 2 of Ref. [IJ] that at the magic 813 nm wave- 
length for Yb, trapping of ^■^'^Cs is magic for B = 2.9G. 

In another approach one may choose to operate on 
multiphoton Mp -^ —Mp transitions [5]; such transitions 
have "magic" values of B fields, where the Zeeman sen- 
sitivity is removed at points where dv/dB = 0. Again 
there is a range of lattice wavelengths where the trap- 
ping may be made magic by varying the angle between 
the lattice wavevector and the quantizing B-field. For ex- 
ample, for ^•^■^Cs |4, 3) — > |3, — 3) transition, the accept- 
able ranges of magic wavelengths in nm are 898-1591, 
863-880, and 512-796 (see Table I of Ref.t5,]). Since the 
magic wavelengths of Ca and Yb fall within these ranges, 
these species can be co-trapped with Cs. 



TABLE I: Values of Al "magic" angles at "magic" wave- 



lengths for divalent atoms. Values for A„ 



from Ref. h7 





Amagic[^??T'J 


9p (linear pol.) 


9k (circ. pol.) 


Hg 


362 


58.89° 


46.95° 


Zn 


416 


48.10° 


70.82° 


Cd 


419 


48.82° 


68.62° 


Mg 


466 


51.91° 


60.73° 


Ca 


739 


52.90° 


58.54° 


Yb 


759 


52.91° 


58.52° 


Sr 


813 


52.93° 


58.47° 



An important consideration is the sensitivity of transi- 
tion frequency to experimental error in the magic angle. 



VI. CONCLUSION 

We have presented a method to simultaneously remove 
Stark shifts for two atomic species trapped in a shared 
optical lattice. One transition in each atom may be made 
insensitive to the intensity of the trapping laser fields. 
This method could prove useful for quantum informa- 
tion and precision measurement experiments. For ex- 
ample, two different atomic species may be used for exe- 
cuting quantum multi-particle gates and state-dependent 
transport (see, e.g. [l3|). Operating the common lattice 
at the simultaneous magic trapping conditions would re- 
move detrimental Stark-induced decoherences for both 
types of qubits. 
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